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ABSTRACT 

Unusual extinction curves of high-redshift QSOs have been taken as evidence that dust is primarily 
produced by supernovae at high redshift. In particular, the 3000 A Todini-Ferrara-Maiolino kink in 
the extinction curve of the z = 6.20 SDSS J1048+4637 has been attributed to supernova dust. Here 
we discuss the challenges in inferring robust extinction curves of high-redshift QSOs and critically 
assess previous claims of detection of supernova dust. In particular, we address the sensitivity to 
the choice of intrinsic QSO spectrum, the need for a long wavelength baseline, and the drawbacks in 
fitting theoretical extinction curves. In a sample of 21 QSOs at z ~ 6 we detect significant ultraviolet 
extinction using existing broad-band optical, near-infrared, and Spitzer photometry. The median 
extinction curve is consistent with a Small Magellanic Cloud curve with Ah^q '^ 0.7 mag and does 
not exhibit any conspicuous (restframe) 2175 A or 3000 A features. For two QSOs, SDSS J1044-0125 
at z = 5.78 and SDSS J1030-I-0524 at z = 6.31, we further present X-shooter spectra covering the 
wavelength range 0.9-2.5 /im. The resulting non-parametric extinction curves do not exhibit the 
3000 A kink. Finally, in a re-analysis of literature spectra of SDSS J1048-I-4637, we do not find 
evidence for a conspicuous kink. We conclude that the existing evidence for a 3000 A feature is weak 
and that the overall dust properties at high and low redshift show no significant differences. This, 
however, does not preclude supernovae from dominating the dust budget at high redshift. 
Keywords: accretion, accretion disks — dust, extinction — galaxies: high-redshift — quasars: general 
— quasars: individual (SDSS J1030+0524, SDSS J1044-0125, SDSS J1048+4637) 



1. INTRODUCTION 

Characterizing the properties of dust in galaxies at 
high redshift is important for our understanding of 
the physics of dust formation and its interplay with 
galaxy evolution. The very large dust masses found in 
some high-redshift QSOs from sub-mm and mm detec- 
tions of their far-inf r ared emission (Priddev et al.l 120031: 
Bertoldi et all [20031: iRobson et ahl 120041: iPriddev et all 



Dweket all [20071: 
2011allbl ld). 



2008[ ) pose a challenge to current models of the evolu- 
tion of massive stars and galaxies. This is because the 
time since Big Bang at z > 6 (< 1 Gyr) is very short 
for large amounts (> 10® M©) of du st to be produced 
by stellar or non-stellar so urces (e.g., 
iMichalow ski et al.l}2010a|[bl : iGall etahl 

MaiolinoeLalJ (1200^ observed the z = 6.20 QSO 
SDSS J1048-I-4637 and inferred an unusual ultraviolet 
(UV) extinction curve, inconsistent with known local 
extinction curves for the Milky Way or the Magellanic 
Clouds, but consistent with predictions of simple theo- 
retical models of dust formation in supernova (SN) e j ect a 
(|Todini fc Ferraral[200ll : IBianchi fc Schneideri[2007h (for 
a variety of extinction curves, see Figure[T]). The intrigu- 
ing suggestion is that SNe are responsible for the rapid 
formation of dust at high redshift, whereas Asymptotic 
Giant Branch stars contribute significantly only at lower 
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redshift, when they ha ve had tim e to ev olve. 

More recently, Galle rani et"aD (|2010L hereafter IGIOI ) 
studied a larger sample of high-redshift QSOs, including 
SDSS J1048-I-4637, and fitted their spectra with a range 
of theoretical S N extinction curves, as suming a QSO tem- 
plate spectrum (jReichard et al.l[2003l ) with an adjustable 
power law slope for the underlying QSO continuum. The 
best-fitting extinction curves were found to be slightly 
more shallow than in the Small Magellanic Cloud (SMC) 
but no stron g evidence for extinct ion curves of the kind 
suggested bv lMaiolino et al.| (|2004|) was found. However, 
obtaining robust results from such an approach is chal- 
lenging because of (i) the narrow rest wavelength range 
considered (~ 0.13-0.3 ^m), (ii) the unknown intrinsic 
(unextinguished) QSO spectrum, (iii) the difficulty in es- 
timating the continuum baseline (QSOs are dominated 
by broad and blended emission lines almost everywhere 
in their UV spectra), and (iv) the possible presence of 
variability, which will give rise to spurious features in the 
spectrum if its constituent data points were not obtained 
simultaneously. 

While the prospects of interpreting unusual extinction 
curves as evidence for SN dust at high redshift are excit- 
ing, it is important to note that detecting significant ex- 
tinction requires sufficiently blue intrinsic spectra of the 
corresponding QSOs. Indeed, previous extinction curves 
reported for high-redshift QSOs have required extremely 
blue QSOs. This would be consistent with the notion 
that only very blue QSOs with reddening can be detected 
in the first place because of the detection limit in the 
SDSS z' band. However, it is a concern that no QSOs 
at high redshift with normal intrinsic colors have been 
found to have extinction, especially in view of the delicate 
procedure of inferring an extinction curve. Moreover, 
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Figure 1. Absolute (bottom) and relative (top) extinction curves. 
Extinction curves are colored and normalized to A1450 = 1 mag 
(this corresponds to Ay = 0.195 mag for the SMC curve): Milky 
Way (blue), SMC (green), Todini-Ferrara theoretical extinction 
curve for SN dust (red). The Todini-Ferrara extinction curve is 
only plotted in the bottom panel because it is not defined above 
3200 A. Artificial extinction curves are shown in gray and are not 
normalized. They are obtained from a thermal disk assumed to be 
a power law for T* = 40, 50, 60, 70, 85, 120, 190 X 10^ K (from top 
to bottom), as discussed in Section |2. II Note the strong similarity 
to the SMC curve. 

very blue QSOs are more likely to have a s iRnificant con- 
tribu tion from a thermal accretion disk (jPerevra et alJ 
|2006[) which is different from the usually assumed intrin- 
sic power-law spectrum. Hence, it would be important 
to test if intrinsic slopes are indeed very blue and, if 
so, to probe fainter QSOs which would allow extinction 
measurements in systems with normal intrinsic spectral 
slopes. 

There is only one example so far of a QSO, SDSS 
J1048-I-4637, with an observed extinction curve ex- 
hibiting a character istic bump at around 3000 A 
(jMaiolino et al ."2004') simil ar to the th eo retical extinc- 
tion curve of Todini & Fc rraral (|2001l ). iStratta et al.l 
11007) did claim evidence for SN dust iii GRB 050904, 
but this was l ater refuted by [Zafar et all (120101), as ac - 
know ledged by 



Stratta et alt (1201 ID. iPerlev et al.l (IMol) 



and I Jang et al.l (|201lD have reported a lTodini fc Ferraral 
(J2001D like extinction curve in GRB 071025, which how- 
ever relies on a photometric redshift and hinges on the 
prec ision of the calibration of the iJ-band measurements. 
IGIOI did interpret the somewhat shallower UV slopes 
found in their extinction curves as being consistent with 
SN dust, howev er, of a kind with no features in the ex- 
tinction curves (jHiras hita et al.l[2008[ ). It remains to be 
proven that an extinction curve with a UV slope, which 
is slightly shallower than the average SMC curve, would 
necessarily require SN dust, considering that radiative 
transfer effects will also give rise to shallower slopes 
and in view of the diversity of U V extinction curves 
in the Milky Way and the SMC (ICardeUi et al.1 [l98l 



iFitzpatrick fc Massal [2007t iGordon et all I200I . Thus, 
observationally, unusual extinction curves are in need of 
confirmation. 

Regarding theoretically predicted ex tinction curves of 
SNe, there are considerable differences (Todini & Ferrard 
20011: iBianchi k Schneider 2007 : Hirashita ct al. 20081 
The underlying dust formation models are sensi- 



tive to assumptions about the parameters controlling the 
dust g rain properties (such as morphology, s ize or compo- 
sition, iGah et al.ll2011cl : iFallest et al1l2011| ) which deter- 
mine the shape and characteristics of extinction curves 
(see Section 16. 5p . Furthermore, dust grains produced 
by SNe will be subject to either disruptive, destructive 
or growth processes due to, for example, shock interac- 
tions in the SN remnant (e.g. , 'Bianchi fc Schneidcrl l2007l: 
Nozawa et~all l2007t iDunne e t al. 2009; Matsuur a et all 
20111) or the interstella r medium (e.g. 



Michalowski et al1l2010al ). or due to reprocessing by the 



intense UV radiation in star-forming regions. Any mod- 
ifications of the original dust grains formed by SNe will 
lead to changes in either the mineralogy or the grain-size 
distribution of the dust, and hence the resulting extinc- 
tion law (e.g., Hirashita & Kuo 2011). It is therefore 
not obvious to what extent the theoretical predictions of 
dust extinction curves from SN models are directly rel- 
evant to extinction curves inferred from lines of sight to 
QSOs. Therefore, it would be prudent not to rely on spe- 
cific parametrizations when inferring extinction curves 
in high-redshift QSOs. A non-parametric approach is 
desirable because we cannot assume that current model 
predictions are correct or directly applicable to observa- 
tional data. 

Thus motivated, this paper is devoted to a study of 
possible extinction signals in high-redshift QSOs ob- 
served by the Infrared Array Camera (IRAC) on the 
Spitzer Space Telescope. We first discuss the procedures 
of inferring QSO extinction, including the choice of in- 
trinsic spectrum, specifically power-law models vs. ac- 
cretion disk models (Section 2). Next we show that the 
long- wavelength IRAC/ Spiteer data allow us to constrain 
the intrinsic spectral slopes of the QSOs and use these to 
determine rough broad-band extinction curves for the full 
sample of QSOs (Section 3). In Section 4 we present X- 
shooter observations of high-redshift QSOs, which, com- 
bined with IRAC/ Spitzer data, allow us to detect ex- 
tinction in high-redshift QSOs with normal colors and 
little extinction. The targets are the broad absorption 
line (BAL) QSO SDSS J1044-0125 at z := 5.78 with a 
large inferred dust mass and the (non-BAL) QSO SDSS 
J1030-f 0524 at z = 6.31 with no sub-mm detection. Dust 
extinction has not been detected in these systems previ- 
ously (GIO). Finally, we study the prototypical system 
for a special SN extinction curve, SDSS J1048+4637, and 
show that there is no strong evidence for the previously 
claimed kink in the extinction curve (Section 5). We 
discuss our results in Section 6. 

2. INFERRING DUST EXTINCTION CURVES IN 

QSOS 

In this section we address the problem of inferring dust 
extinction in QSOs from the wavelength-dependent at- 
tenuation of the emitted light. The standard method 
of obtaining an extinction curve is to measure the ratio 
between an extinguished spectrum and that of a refer- 
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ence spectrum, typically an unextinguished object or a 
model spectrum. For example, multiply imaged (gravi- 
tationally lensed) QSOs have been used successfully to 
measure extinction curves o f intervening (lensing) galax- 
ies (Jaunscn & Hiorth 1997; Falco et al. 1999: 'Toft et alJ 
[2000 : Motta ct al. 2002 : Eliasdottir ct al. 2006)) and 
power-law forms to measure ex tinction curves u sing 
gamma-ray burst afterglows ('e.g.. lZafar et al.|[2010[ ). 

However, to measure extinction curves from the spec- 
tra of unlensed QSOs one must resort to comparing the 
observed spectrum with a model intrinsic spectrum. The 
interpretation of any signal may strongly hinge on the as- 
sumed shape of the intrinsic spectrum. 

It is believed that a thermal accretion disk is respon- 
sible for the majority of the emission in the UV-o ptical 
regime of QSOs ' (Shang et al. 2005; Pcrc vra et al.l (20061 : 
iKishimoto et al.l 120081 : ILawrencd 120121 ). On the other 
hand, it is standard practise to describe the contin- 
uum emission in QSOs as a power law, /i, ex i^", i.e., 
f\ oc A~^~", presumably resulting from reprocessing of 
the flux from the thermal disk. 

A thermal disk model is defined as 
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(IShakura ~Sunvaevl 11971 IShieldsl [1971 IFrank et al.l 
[l992i : iPerevra et al.l 120061 ). In the above expressions, G 
is the gravitational constant, c is the speed of light, Mbh 
is the black-hole mass, h is Planck's constant, k is Boltz- 
mann's constant, A/acc is the mass accretion rate, crs is 
the Stefan-Boltzmann constant, and s is the radius di- 
vided by the innermost disk radius. The latter is taken 
to be the last stable orbit for a non-rotating black hole, 
i.e., 3 times the Schwarzschild radius. The characteris- 
tic temperature T* uniquely controls the shape of the 
thermal disk spectrum. For AT* — > oo the spectrum 
asymptotically approaches an a = 1/3 power law. 

In Figure [2] we show a thermal disk model as well as a 
power-law model for the continuum of a template QSO 
spectrum. It is evident that a power-law model (with 
a = —0.46) provides a better description of the underly- 
ing QSO spectrum (i.e., unaffected by emission lines) at 
wavelengths redward of Lya. A thermal disk model with 
T* ss 85000 K provides a reasonable approximation, but 
overshoots the continuum of the template spectrum in 
the range 0.13-0.23 ^m. To avoid overshooting requires 
T* < 60000 K, but such models severely underrepresent 
the UV portion of the spectrum. 

2.1. Artificial extinction curves from accretion disk 
spectra modeled as power laws 

If the assumed intrinsic model spectrum is wrong, then 
the method of obtaining extinction laws relative to this 
model will lead to spurious signals. For example, if the 
intrinsic spectrum is not a power law but modeled as 
such, results mimicking standard extinction laws could 
be inferred, even in the absence of dust. To quantify 



100 



CO 




Wavelength ()j,m) 



Figure 2. QSO template spectrum and models. Standard ther- 
mal disks with a range of characteristic temperatures are plotted as 
gray curves (T* = 40, 50, 60, 70, 85, 120, 190 X 10^ K from bottom 
to top, solid curves and T* — > oo, dashed line). Also plotted is an 
H5" r/SDSS template spectrurn (violet, based on Tclfcr ct al. (200^) 
and I Vanden Berk et al.l 11200 ID ). None of the thermal disk models 
provide acceptable representations of the template spectrum con- 
tinuum. For example, the T* ^ 85000 K model overshoots the 
continuum of the template spectrum in the range 0.13-0.23 fira. 
Between 0.13 fim and 0.50 fim a power law with spectral index 
a = —0.46 (green line) is a better fit than the thermal dis k, con- 
sider ing the effects of the broad emission lines ( Vanden Berk et al.l 
120011) . The SDSS template spectrum is not plotted above the loca- 
tion of [O III] A5007 where it is likely contaminated by host galaxy 
emission. This is not expected to be an issue for the high-redshift 
QSOs studied in this paper. 



the effect of a thermal contribution to the UV-optical 
spectrum of a QSO we compute the artificial extinction 
signal inferred when assuming a power-law intrinsic spec- 
tral distribution for a continuum described by a standard 
thermal accretion disk. 

To obtain an effective extinction curve from a thermal 
model we (wrongly) assume the thermal accretion disk to 
be a power law between 5150 A and 8300 A (to mimic the 
analysis in the following sections) and compute the ratio 
between the accretion disk spectrum and the power law. 
The results are shown in Figure [TJ We retrieve featureless 
'extinction laws' which are similar to, but not identical 
to, the SMC extinction curve. The inferred 'extinction' 
is anti-correlated with the intrinsic temperature of the 
disk. This effect must be kept in mind when interpreting 
observational results although we recall that Figure [2] 
demonstrates that the contribution from such a thermal 
spectrum cannot be dominant. 

For the analysis of observational data, we will therefore 
assume a power-law continuum model between Lya and 
up to 8300 A. This is not a good assumption at longer 
wavelengths where hot dust emission has been shown 
to dominate th e emission in most high-redshift QSOs 
(| Jiang et al.ll20 6, 20 10). We discuss hot dust contami- 
nation in Section F3.2l and quantify the effect of a thermal 
disk contribution on our extinction measurements in Sec- 
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Figure 3. Optical/NIR and IRAC /Sp i tzer b road-band data for 
21 ^ ~ 6 QSOs from Jian g et al.1 ll200g . [20Tct) . Tiie data points 
are not plotted at fixed rest wavelengtlis because of tiie range in 
redshift (5.78 < z < 6.42). Note the conspicuous effects of tiic 
emission lines of Ly« and Ha in the z and 4.5 fim bands, respec- 
tively. For clarity, error bars are not plotted. The typical photo- 
metric uncertainties are 10-20 % in the optical/NIR bands and for 
the 5.8 and 8.0 fira IRAC bands and at the percent level for the 
3.6 and 4.5 ^m data points. 



tion 13.31 

3. Spitzer SAMPLE OF Z ~ 6 QSOS 

We used IRAC /Spitzer photometry of 21 z ~ 6 QSOs 
supplemented by tabulated (but non-simultaneous) op- 
tical and near-infrared bro a d-ban d photometry pub- 
lished by iJiang et all ()2006L l2010t ). The optical and 
near- infrared d ata were corrected for Galactic extinc- 
tion (jSchlegel e t al. 1998) assuming Ry = 3.1. These 
corrections are small, with E{B — V) ranging from 0.01 
to 0.1 mag with the exception of SDSS J0353+0104 for 
which E{B — V) — 0.29 mag. The median value is 
E{B — V) = 0.03 mag. Such small Galactic extinctions 
are reassuring because the possible systematic biases on 
the QSO extinction signals from foreground extinction 
removal are thereby minimized. 

In Figure [3] we show the resulting broad-band spectral 
energy distributions. We do not plot the 24 /xm data 
points because they are dominated by hot dust emis- 
sion in most cases. The 3.6 ^m, 4.5 /xm, and 5.8 iim 
data points do not seem to be affected (jJiang et al.ll2006l 
|2010| ). The 8.0 /J,m data points are in some cases con- 
taminated by hot dust emission. For this reason, these 
data points are not included in our analysis but shown 
in Figure [3] to illustrate the effect. Unfortunately, the 
4.5 ^m photometry is significantly affected by Ha emis- 
sion at z ~ 6 and the 3.6 /im less so by H/3 and [O III] 
AA4959, 5007 emission (see Figure[2l). 

3.1. Correction for emission-line contamination 

We correct the IRAC photometry for contamination 
by broad emission lines as follows. It is evident from 
Figure [3] that while there is a clear bump in the 4.5 
/zm band there is no conspicuous effect in the 3.6 /im 
band. For each QSO we therefore fit a power law to 
the 3.6 nva and 5.8 /xm fluxes and estimate the excess 



Figure 4. Relation between Ho to H/3-|-[0 III] equivalent width 
ratio and luminosity of the QS Os. The da ta points arc for z < 
0.4 Palomar-Green QSOs (Shang et al.''2007). The SDSS QSO 
template ratio (Vanden Berk et al. 2001) is indicated as a dashed 
line spanning the range of SDSS QSO absolute magnitudes. The 
dotted line indicates the ratio adopted in our analysis, spanning 
the range of absolute magnitudes of the z ~ 6 QSOs. 

Ha flux. The median excess flux is 25%. correspond- 
ing to an Ha equivalent width of ^ 460 A. This value 
is larger than for t he template SDSS spectrum (194 A; 
iVanden Berk et an r2001') and consistent with the median 
value for the complete samp le of z < 0.4 Palomar-Green 
QSOs (434 A; Shang et al.ii2007[ ). This strongly suggests 
that the bump is indeed due to Ha emission and sup- 
ports our method of analysis. We then correct the 3.6 
/xm band using an assumed equivalent width emission- 
line contamination rat io, 7? = Ha/(H/3-|-[0 III]). The tem- 
plate SDSS spectrum (jVanden Berk et aLll2001[ ) suggests 
77 = 3.5 while the brighter Palomar-Green QSOs yield a 
median 77 = 4.2, ranging from 2.7 to 5.3 (|Shang et al.l 
l2007f l. with a trend for more luminous QSOs to exhibit 
larger 77 (Figure 2]) . These QSOs are on average an or- 
der of magnitude brighter than the z < 0.4 SDSS QSOs 
making up the region around Ha and are therefore likely 
more relevant for our purposes. According to Figure SI 
a ratio of 77 = 4.75 appears to be a reasonable value for 
luminous QSOs, which we adopt in what follows. We 
note however that an even higher ratio could in principle 
be relevant for some of the highly luminous z '^ 6 QSOs. 
Taking into account the relative wavelength ranges of 
the relevant IRAC filters gives a ratio of 3.5 between the 
Ha to 11/3+ [O III] emission-line to continuum contamina- 
tion, which we apply to the 3.6 ijm flux. The procedure 
was repeated by refitting a power law to the corrected 
3.6 /im and 5.8 /im fluxes until convergence and resulted 
in a median correction of —8% to the 3.6 /im flux. The 
uncertainty in the correction induces a small systematic 
uncertainty in the inferred intrinsic power-law slope, a, 
which we address below. 

3.2. Extinction 

The relative extinction curve, ^(A), for each QSO is 
obtained as the ratio between the restframe UV photo- 
metric data points and the extrapolated power-law model 
restframe optical spectrum, corrected for emission lines 
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Figure 5. Broad-band extinction curves relative to the 

IRAC/ Spitzer power-law baseline for 21 z ~ 6 QSOs. The thick 
black curve, with an indicative error of 0.20/\/21 = 0.04 mag, is 
the median extinction curve of the 21 QSOs while the gray curve 
is a representative SMC extinction curve with A1450 = 0.33 mag. 
Note that for the SMC, Ay ~ O.4IA1450. A 2175 A bump, charac- 
teristic of the Milky Way extinction curve (dashed gray curve), is 
not seen in the median extinction curve. There is no conspicuous 
3000 A kink either. 



as explained above. Figure [5] shows the resulting ex- 
tinction curves. There is significant scatter in the data, 
which are affected by other emission lines, depending to 
some extent on the actual redshift of each QSO. Vari- 
ability may also play a role in some cases. Finally, there 
are substantial photometric uncertainties associated with 
the optical and near-i nfrared photometr y corresponding 
to about 0.1-0.2 mag (jJiang et al.ll2010f ). 

Nevertheless, it is evident that the curves do not scat- 
ter around zero. If the data were dominated by intrinsic 
variability, they would scatter symmetrically around the 
zero mean. This is not the case. There is a clear positive 
effect, as expected if the QSOs are reddened relative to a 
power-law spectrum. The effect is relatively small, how- 
ever, with values reaching a maximum of about 1 mag 
in the restframe UV. Only one out of 21 QSOs (SDSS 
J0840-I-5624 at z = 5.85) appears to exhibit a significant 
negative signal. This may be due to variability between 
the epochs at which the SDSS photometry, the near- 
infrared photometry, and the IRAC / Spitzer photometry 
were obtained. Alternatively, since we are relying on de- 
riving a slope from just two literature data points, we 
cannot exclude problems with the reported photometry. 
In fact, the 8.0 /im data point for this target seems to 
indicate a somewhat steeper spectral slope which would 
give rise to a less negative or even positive extinction 
signal. 

The data are not of sufficient quality to obtain the 
detailed shapes of the individual extinction curves. The 
median extinction curve is consistent with an SMC curve 
with Av = 0.13 mag while a Milky Way 2175 A bump 
is strongly ruled out. Likewise, there is no conspicuous 
3000 A kink. 

For each QSO we can obtain an estimate of the UV 
extinction, A1450, from the 1450 A magnitude (the z and 
J bands are affected by Lya and C IV emission, respec- 
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Figure 6. Spectral slopes versus extinction (A1450, violet sym- 
bols) for 21 IRAC/Spitzer 2 ~ 6 QSOs. Also shown are the re- 
sults of jGlQl (green sym bols, stars for detections, arrows for up- 
per limits) . We used the IGIOI appropriate mean extinction curve 
or the extinction curves of individual objects to obtain j4i45o 
from their A3000 a-^d assumed the SMC extinction law relation 
^1450 = 0.47yli45o to obtain A1450. Our corresponding results 
for the three green stars are shown as violet stars. The dashed 
and dotted lines indicate the degeneracies between A1450 and a 
discussed in the text. The yellow region represents an exclusion 
zone due to the selection effect that reddened QSOs with a given 
IRAC 5.8 /im flux are more easily detected when the QSO is blue. 
The borderline between the excluded (yellow) and allowed region 
represents a QSO with an intrinsic a = spectrum, detectable 
at 10 cr at 5.8 fim and observable to z' = 20.9 mag in the SDSS. 
The gray line indicates the relation between excess inferred red- 
dening and inferred spectral slope for a pure thermal disk, ranging 
from values around a Ri — 0.2,y4i45o ~ 1.3 for low temperatures 
(T* Si 40000 K) to the asymptotic values (a = 0.33, A1450 = 0) at 
high temperatures. 

tively) . Figure [S] shows the relation between the inferred 
extinction, Ah^q, and spectral slope, a. There is a clear 
correlation, which can be understood as follows: For a 
power, law f^ o: v", anchored at Aanchor and observed at 
Aobs the difference in the inferred relative extinction at 

Aobs IS 

5A^2.5log{^^^^^^^]Sa (3) 



^anchor 
Aobs 

if the spectral slope changes by 5a. 

One uncertainty in a arises from the slightly model 
dependent correction for the flux of H/3 and [O III] in the 
3.6 fiia band. The spectral slope is anchored at the 5.8 
/im IRAC band, i.e., at around 8300 A in the restframe 
at 2; = 6. The observed wavelength is at 1450 A in the 
restframe. Hence, we expect the inferred extinction and 
spectral slope data points to be degenerate along 

6 A = l.89Sa, 

plotted as a dashed line in Figure El 
variation Srj leads to a change in a of 

4.75' 
which indicates that this is a small effect. 

Another source of error is the possible contamination of 
the 5.8 fim band by hot dust emission. If this uncertainty 
is dominant, the spectral slope is effectively anchored at 
the 3.6 nui IRAC band, i.e., at around 5150 A in the 



(4) 
We find that a 

(5) 
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restframe at z = 6, and 



5A = 1.375a. 



(6) 



This relation is plotted in Figure [Has a dotted line. We 
find that the median value of a in this case varies as 

Sa^-l.S^-^^^i^, (7) 

which shows a strong sensitivity to the 5.8 /Km flux. For- 
tunately, the hot dust contamination at 8.0 /xm is small 
(< 10%, see Figure El and hence negligible in most cases 
at 5.8 /im: the Planck function drops by a factor of 9 
between 8.0 /im and 5.8 /im in the Wien region, while 
the flux is roughly constant for a k,Q. This is consistent 
with the hot dust models of IJiang et al.l (J2006D which in- 
dicate that the contamination of the 5.8 /xm flux is of 
order 1%. We note that any contamination by hot dust 
emission will cause_the 5.8 /xm flux to be overestimated 
and hence a and A to be underestimated. For consis- 
tency, we applied a —1% correction to all 5.8 /xm fluxes 
to account for contamination by hot dust emission. 

Having shown that the systematic uncertainties related 
to the choice of rj or due to hot dust contamination are 
small, we conclude that the actual photometric uncer- 
tainty in the 5.8 /im flux dominates the error budget 
(the uncertainties in the 3.6 and 4.5 /Ltm fluxes are much 
smaller). The median relative error in the 5.8 /im flux is 
8% which translates into a typical error in A1450 of 0.20 
mag, or 0.04 mag in the mean for the sample of 21 QSOs. 

Selection effects are described by relations similar 
to Equations ([4]) and (jO]). As the selection of high- 
redshift QSOs is based on a z'-filter detection and i'- 
filter dropout, the dominating selection effect is related 
to the z'-filter brightness. For a given IRAC 5.8 /xm band 
detection, the corresponding constant of proportionality 
in Equation ^ becomes slightly larger, around 2 (tak- 
ing the z'-band to be the observed wavelength), i.e., it 
defines an exclusion zone above a line with a slope of 
5^i45o/'5a = 2 where QSOs are too faint for detection. 
The normalization depends on the intrinsic brightness of 
the QSOs and the depth of the IJiAC / Spitzer observa- 
tion. In Figure[n]we illustrate the effect for a typical QSO 
detected at 100 ± 10 /xJy in the 5.8 /xm band. Assum- 
ing a fiat spectrum (a = 0) this corresponds to a QSO 
with an unextinguished brightness of z' = 18.9 mag. The 
faintest QSO in the sample has z' = 20.9 mag, i.e., we 
can observe_2 mag of relative extinction for a QSO with 
a = 0. iGld discussed a color selection effect based on 
z' — J but this constraint is weaker than the brightness 
constraint in our case and is not plotted. 

3.3. Discussion 0/ the broad-band results 

Figure |6] illustrates the fact that available detection 
limits do not allow the detection of very reddened QSOs 
with small spectral slopes (i.e ., int rinsically red QSOs). 
Therefore, as pointed out by IGIOI . the fact that previ- 
ously reported high-redshift QSOs with significant ex- 
tinction are intrinsically very blue appears consistent 
with the selection effects. On the othe r han d, our re- 
sults appear to be inconsistent with the iGld results. In 
Figure [7] we compare the spectral slopes obtained for the 
QSOs studied here with those obtained bv IGld . There 
is no obvious correlation between the respective values 
and no systematic bias either. It seems the outcomes 
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Figure 7. Comparison of inferred spectral slopes for 2 ~ 6 QSOs. 
We plot the intrins ic sp ectral slopes obtained in this work and 
comp are to those of lGlOl and those adopted bv IJiang et al.l I I2006I , 
120101) . The filled stars represent QSOs for which a is reported by 
the different sources. The symbols should lie on the solid lines 
for consistency between the different methods for inferring the in- 
trinsic spectral slopes. Note that our slopes and those of Gffl are 
extinction corrected whereas those of lJiang et al.i (|2006| , |2010| ) are 
not. This partly explains why the latter values are smaller, corre- 
sponding to redder sources. 



of the two approaches are very different. Hence, given 
the extinction-spectral slope correlation, extinctions ob- 
tained in at least one of the approaches may be unre- 
hable. In the s ample of 3814 SDSS QSOs studied by 
iReichard et al.l (|2003D . QSOs with observed slopes o f 
a ~ 1 or bluer are extremely rare. iDavis et al.l (|2007D 
find an even more restricted range of slopes in the 1450- 
2200 (2200-4000) A range from 3646 (2706) quasars with 
redshifts in the z = 1.67-2.07 (0.76-1.26) range. QSOs 
with a > .5 are practically non-existent in this sample. 
Hence, for IGld to find three extremely blue (a > 0.6) 
z ^ 6 QSOs seems unlikely. None of our values of a are 
this extreme. Our median value of a « —0.2 may be 
compared to the SDSS template slope of —0.46 and the 
mean slopes of - 0.59 (1450-2200 A ) and -0.37 (2200- 
4000 A) found bv lDaviset all ()2007D . This may indicate 
a slight bias towards blue QSOs, which may partly be 
due to a selection effect because even little extinction 
would make intrinsically red QSOs undetectable. 

As discus sed ab ove, the main weakness of the approach 
adopted by iGld is that it relics on an assumed extinc- 
tion curve, fitted over a limited range in wavelength, to 
regions of the spectrum which may not quite reach the 
true continuum level. Some of these weaknesses are elim- 
inated in our approach, however, at the expense of not 
determining absolute extinction values: We emphasize 
that the extinction curves obtained here are relative to 
the 3.6-5.8 /im power law model and in this sense are not 
absolute but strjctly lower limits (for an SMC extension 
above 3000 A, ili45o w 0.5^1450, see Figure [S]). 

As such, our approach is insensitive to a possible 
gray (i.e., wavelength independent) extinction compo- 
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nent. However, for an entirely gray extinction curve our 
approach would have led to zero relative extinction to be 
detected. The fact that there is a clear positive detection 
implies that the typical extinction curve of high-redshift 
QSOs is wavelength dependent, giving rise to reddening 
as expected (still under the assumption of an intrinsic 
power- law spectrum). 

To address the question of systematic errors in the de- 
rived extinction values due to a contribution of a thermal 
disk we overplot the extinction-slope relation inferred for 
a non-reddened thermal accretion disk model in Figure[Bl 
The curve reflects a range of characteristic temperatures 
starting at T* = 40000 K (a = -0.2,ii45o = 1-3 mag) 
and approaching the asymptotic value (a = 0.33, A1450 = 
0) for large T*. Remarkably, this curve is practically 
perpendicular to the degeneracy relations. Equations Q 
and ([5]), so basically adds to the scatter in this relation. 
Moreover, it is evident that the results do not span the 
full range of temperatures plotted, preferring values in 
the range T* = 60000-100000 K. Interestingly, the fact 
that the majority of the points do not lie on this relation 
shows that the QSO spectra are not well described by 
pure thermal disks and require a significant additional 
component, such as the power law assumed in our anal- 
ysis. 

4. X-SHOOTER OBSERVATIONS 

4.1. Targets 

4.1.1. SDSS J1030+0524 

SDSS J103 0-H0524 (IFan et al.ll200l is a z = 6 309 ± 
0.009 QSO (IJiang et al l 120071 iLro muro et al.l [200I . 
iPriddev et al.l (|2003D obtained an upper limit to the dust 
mass of 1.4 X 10^ M© from SCUBA observations. The 
source is located in a crowded sub-mm field. The sub- 
mm sources in th e field are, however, likely due to fore- 
ground galaxies (jPriddev et al.l I2008D . A deep X-ray 
spectrum did not reve al evidence for X-r ay absorption 
Nh < 8 X 10^2 cm-^: llWah et al.|[200l . IJiang et al.l 
2006[) find a hot dust mass larger than 17 Mq. There 
are Lya emitters in the field at z = 5.7, one consistent 
with being a foreground absorber (jGonzalo Diaz et all 
120111) in the QSO spect rum. There is also a s trong Mg II 
absorber at 2; = 2.780 (jSimcoe et al.ll20lll) . IJiang et all 
(|201Q) report a black hole mass of 2.75 x 10^ Mq. 



4.1.2. SDSS JIO44 
J1044-0125 (IFan et al.l 
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10** Mq (IPriddev et al.l I2003D an d a hot 
dust mass larger than 53 Mq (jjiang et al.ll2006D (second 
highest in the sample). It could be magnified by a t most 
a fac t or of 2 due t o gravitational lensing (Shiova et al.l 
I2002t iWvithd I2004D so the d erived prope r ties ar e likely 
not strongly overestimated. IJiang et "alj (|2010D find a 
black hole mass of 10.5 x 10^ Mq. 

4.2. Observations 

We observed SDSS J1044-0125 and SDSS J1030-t-0524 
on 22 and 23 December 2009 UT, respect ively , with 
the X-shooter spectrograph (jD'Odorico ct al.l 120061 : 



iVernet et al.ll20l"ll) mounted at the Cassegrain focus of 
the Kueyen unit of the Very Large Telescope (VLT) 
at the European Southern Observatory (ESO), Cerro 
Paranal, Chile. The observations started at around 7 UT 
on both nights, and consisted of a set of 4 x 1200 s integra- 
tions, with the objects being offset along the slit between 
the exposures. The seeing as measured from the guide 
stars used for the QSO observations was approximately 
0.7" for both nights. During the SDSS J1044-0125 ob- 
servations thin clouds were scattered across the sky, while 
the sky was clear during the next night when SDSS 
J1030-f0524 was observed. The approximate airmass 
range for both QSO observations was 1.1-1.3, i.e., with 
a maximum zenith angle of 40 degrees. 

X-shooter allows for simultaneous spectroscopic obser- 
vations at intermediate resolution in three different arms: 
the ultra-violet and blue (UVB), visual (VIS) and near- 
infrared (NIR) wavebands, with a continuous spectral 
coverage of 0.3-2.5 /J,m. In this paper we only report 
results for the VIS and NIR arms as the UVB arm cov- 
ers restframe wavelengths much smaller than Lya (1216 
A). We used slit widths of 0.9" for both the VIS and 
NIR arms, resulting in resolving powers of 8800 (VIS) 
and 5600 (NIR). Immediately following the QSO inte- 
grations, a telluric standard star was observed at a very 
similar airmass and with the same instrument set-up as 
the QSOs, i.e., with the same slit widths and binning (no 
binning in the spatial direction, and two pixels in the dis- 
persion direction for the VIS arm). During both nights, 
observations of spectrophotometric standards were also 
secured. The QSOs and the standards were all observed 
with the slit aligned with the parallactic angle. 

Using version 1.2.2 of the X-shooter pipeline (see 
iD'Odorico et al.l 1200 6^ in polynomial mode, the object, 
telluric and standard star spectra were reduced to two- 
dimensional bias-subtracted, fiat-field corrected, order 
rectified and wavelength calibrated spectra in counts. 
One-dime nsional spect ra were optimally extracted from 
these fsee lHornill986l) with a custom IDL routine. 

We used the general ve rsion of the I DL-b ased xtell- 
cor program developed by iVacca et al.l ([2003) for SpeX 
- a 0.8-5.5 /zm medium-resolution cross-dispersed spec- 
trograph at the NASA Infrared Telescope Facility - to 
perform a telluric correction and flux calibration of the 
QSO spectra. Xtellcor uses a high-resolution model spec- 
trum of Vega that is scaled, reddened and convolved with 
the instrumental resolution to match the observed spec- 
trum of a telluric standard as close as possible, including 
the hydrogen absorption lines. Comparison between the 
model and observed telluric spectra then allows for cor- 
rection of the telluric features in the object spectra, and 
also provides a flux calibration if the B and V magni- 
tudes of the telluric standards are known. The telluric 
spectra that we used were both observed immediately fol- 
lowing the QSO: HIP 49704 for SDSS J1044-0125 and 
HIP 055011 for SDSS J1030-H0524, at airmasses of 1.17 
and 1.11, respectively; this is within 0.1 airmass of that 
of the QSO spectra. Both telluric standards are clas- 
siflcd as stellar type B9, and are therefore very similar 
in type to an AOV star that is optimally suited for this 
telluric correction method. The advantage of using the 
telluric spectrum to also perform the flux calibration - 
using the Vega model - is that the telluric standard was 
taken immediately following the QSO, at a similar air- 



Hjorth et al. 



mass, and with the same sht width, which minimizes sht 
loss differences. 

We checked the flux cahbration by applying the sensi- 
tivity function (which converts measured counts to ab- 
solute flux) inferred by the xtellcor program directly to 
a spectro-photometric standard star observed during the 
same night, and by comparing the resulting spectrum to 
tabulated flux values. For the SDSS J1044-0125 night 
we used observations of BD-l-17 4708 at airmass 2.63, 
while for the SDSS J1030+0524 night we used GD71 at 
airmass 1.45. These standards were observed with a 5" 
wide slit. On both nights, we find that the inferred flux 
for these standards is a factor of 1.3 higher than the 
tabulated values, while this factor is constant across the 
spectra to within 5% - the approximate accuracy of the 
relative flux calibration. The offset of 30% could partly 
be attributed to slit losses (the telluric standards and 
QSOs were taken with a 0.9" slit). 

Given these absolute flux calibration offsets, we com- 
pared the resulting flux-calibrated spectra for the QSOs 
with the JHK m agnitudes collected in Table 1 of 
I.Tiang et all (|200l). For SDSS J1044-0125 the spec- 
trum flux values are consistent with the magnitudes 
within their errors (0.05-0.1 mag). However, for SDSS 
J1030-f0524 the spectrum flux is too high by 30% as 
measured from the H and K bands where the f\ is flat 
and void of spectral features. 

This difference prompted us to check if the flux of ei- 
ther QSO is variable in time. We measured the magni- 
tude of SDSS J1030-t-0524 from X-shooter z'-band acqui- 
sition images using aperture photometry and calibrated 
against the SDSS. We find z' = 19.81 ±0.06 on the night 
of 23 December 2009, while we find z' = 19.93 ± 0.04 
for the night of 8 February 2011 (taken from the ESO 
archive), i.e., reasonably consistent with being constant. 
This compares to z' = 20.05 ±0 .10 as reported in the dis- 
covery paper of|^eiin (gOOil) . For SDSS J1044-0125, 
we measure a magnitude of z' = 19.2 7 ± 0. 07 from 
our acquisition images, while iFan et aTl (|2OOO0 report 
z' = 19.23 ± 0.07. No other X-shooter acquisition im- 
ages, or a suitable series of images taken with other ESO 
instruments are available to check for flux variability of 
SDSS J1044-0125. The broad-band photometry is con- 
sistent with a marginally significant brightening of SDSS 
J1030+0524 of 0.24 ± o'.12 mag between 2000 and 2009, 
while SDSS J1044— 0125 is consistent with being con- 
stant. These findings refiect the inconclusive spectro- 
scopic results. Given that we do not have information of 
the brightness of SDSS J1030+0524 at the time of the 
TRAC / Spitzer observations we proceed by not making 
any corrections due to variability but note that a 30% 
brightness correction would offset the resulting extinc- 
tion curve higher by about 0.30 mag. 

Finally, we corrected the fiux-calibrated QSO spec- 
tra for the Galactic extinction along the sight lines, with 
E{B - V) of 0.024 for S DSS J1030±05 24 and 0.051 for 
SDSS J1044-0125 using ISchleeel et all ll998). We dis- 
play the X-shooter/ VLT spectra in Figure |51 

4.3. Extinction curves 

QSO spectra are strongly affected by emission lines, 
both bright conspicuous broad emission lines and fainter, 
blended lines giving rise to a quasi-continuum, especially 
around 0.2-0.3 /zm, the so-called small blue bump due 
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Figure 8. X-shooter/VLT spectra of SDSS J1030+0524 and 
SDSS J1044— 0125. The gaps in the spectra are due to the main 
atmospheric absorption bands between the J and H band and the 
H and K band. Regions plotted in gray are not included in the 
analysis due to the presence of strong emission lines (see Figure |9]l, 
broad absorption, or intergalactic absorption. The main emission 
lines are indicated. The spectra have been smoothed to 0.02 fim 
resolution. 
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Figure 9. QSO emission-line to continuum ratio obtained as the 
ratio between the QSO template and the continuum model shown 
in Figure [2] The solid line represents identity between the QSO 
template and the continuum. Note that the continuum is not 
reached at any point above 1450 A. The dashed line indicates re- 
gions above which the spectrum is severely affected by strong broad 
emission lines. 



to B aimer and Fe II lines. It is essential to remove the 
contamination by these lines to determine the true con- 
tinuum level, ass umed to be a powe r law i n our analy- 
sis. We use the IVanden Berk et al.l (l200ll) SDSS QSO 
template to estimate the ratio between the QSO signal, 
affected by emission lines, and the continuum level (Fig- 
ure [9]). We divided the X-shooter spectra by this ratio 
to obtain an estimate of the quasi continuum, shown in 
Figure [TUl It is evident that our procedure of removing 
the emission line contribution works reasonably well in 
these systems. 
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Figure 10. X-shooter/VLT spectra of SDSS J1030+0524 and 
SDSS J1044— 0125 after correction for emission lines, see Figure|9] 
The filled circles are the IRAC / Spitzer data points. The dashed 
line is the adopted intrinsic power-law spectrum fitted to the IRAC 
data points as described in the text. The extinction is inferred as 
the ratio bet ween the dashed line and the spectrum. As shown 
in Figure [TT1 the slight mismatch between the extrapolated power 
laws and the observed spectra are eliminated when choosing the 
1(7 lower boundary for the 5.8 /im fluxes. 

We repeated the procedure of the previous section 
to tie the spectra of SDSS J1030+0524 and SDSS 
J1044-0125 to the IRAC/Spitzer data. We assume that 
there has been no variability of the QSOs and that the 
intrinsic continuum spectrum is well represented by a 
power law as fixed by the 3.6 /xm and 5.8 /zm IRAC data 
points, suitably corrected for emission lines. We can then 
determine the extinction law as the ratio between the 
extrapolated power law and the emission-line corrected 
X-shooter spectrum. The resulting extinction curves are 
presented in Figure [11] 

The signals are not very strong, although_there seems 
to be the expected reddening effect, with ^1450 = 0.14 
mag for SDSS J1044-0125 and Amso = 0.24 mag for 
SDSS J1030+0524, consistent within the photometric er- 
rors with the broad-band values for these QSOs (Figures 
[5] and [5]). Neither QSO had evidence for extinction in 
the iGlOi analysis, but their limits are fully consistent 
with our detections. The systematic errors in these ex- 
tinction curves amount to ± '-^ 0.25 mag at 1450 A for 
these systems because of the ~ 11% errors in their 5.8 
fiia fluxes. Nevertheless, we can rule out conspicuous 
features at 2175 A (and to some extent 3000 A). SDSS 
J1030-I-0524 has a UV rise which is slightly shallower 
than the SMC, while there is an apparent downturn in 
SDSS J1044-0125. 

5. SDSS J1048+4637 

An apparent UV-flattening was found in the extinction 
curve of SDSS J1048-I-4637, based on prism spectropho- 
tometry using the Near Infrared Camera Spectrometer 
(NICS) instrument at the Telescopio N azionale Galileo 
(TNG ) and attributed to dust from SNe (iMaiolino et al.l 
I2004D . In a reanalysis of the same data, iGld found that 
the presence of such a steep-flat-steep ('kinked') extinc- 
tion curve was heavily dependent on the intrinsic slope 
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Figure 11. Inferred X-shooter/VLT extinction curves for SDSS 
J1030-I-0524 and SDSS J1044-0125. The thick black line is a 
smooth 3rd order polynomial fit to the data. The thick dashed 
lines corresponds to assuming the lower la error on the 5.8 /tm flux. 
If a 30% brightness correction is applied to the SDSS J1030+0524 
spectrum (see Section 14.21 1 the resulting extinction curve would be 
offset higher by about 0.30 mag. Indicative relative SMC (solid) 
and MW (dashed) extinction curves are overplotted in gray. 

assumed. A different intrinsic slope resulted in a rela- 
tively flat extinction curve without a kink. We have at- 
tempted the construction of extinction curves for SDSS 
J1048+4637 using data from the literature, acquired us- 
ing the SAO/NASA Astro physics Data Sys tem's Dex- 
ter data extraction applet (jDemleitner et al][2 001). Op- 
tical/N I R spectra wer e acqu ired from Maiolino et sJ] 
(pool . Ilwamuro et all (|200l . iFan et al . (2003, 200^^ 
based on data from NICS/TNG, from the OH-Airglow 
Suppressor/Cooled Infrared Spectrograph and Camera 
for OHS (OHS/CISCO) on Subaru, from the Low- 
Resolution Spectrograph (LRS) on the Hobby-Eberly 
Telescope (HET) , and from the Echelle Spectrograph and 
Imager (ESI) on Keck. 

These data were all analysed with the method de- 
scribed in Section 14.31 The resulting extinction curves 
using the TNG and the Subaru-(-Keck data, respectively, 
are shown in Figure [12] (we do not show the LRS/HET 
data which are consistent with the ESI/Keck data, but 
of lower quality) . 

The data from NICS/TNG r esult in an ex t inctio n 
curve similar to that found by IMaiolino et al.l ()2004l) . 
There is however no conspicuous bump in the extinction 
curve, although a plateau is evident. In the combined 
ESI/Keck+OHS/CISCO/Subaru spectrum the overah 
normalization of the extinction curve is significantly 
lower and instead of a plateau there is a monotonically 
rising extinction towards the red. The shapes of the ex- 
tinction curves do not seem very different though, with 
only a mild change of slope in the observed NIR. In both 
cases there is a UV upturn with a slope similar to the 
SMC. 

The differences in the derived extinction curves may 
be due to variability and/or there may be instrumen- 
tal causes. In case of variability the source would have 
had to vary by as much as 60% in the restframe UV to 
explain the difference in the absolute extinction level. In- 
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Figure 12. Relative extinction curves for SDSS J 1048+4637. The 
green curve is based on our re-analysis of the NICS/TNG data 
jMaiolino ct al."2004') whereas the violet curve is based on a com- 
bination of ESI/Kcck and OHS/CISCO/Subaru data. The black 
overplottcd curves are 3rd order polynomial fits to the data, for 
illustration. For comparison with the ESI+OHS curve, the rela- 
tive SMC extinction curve is overplotted in gray. The curves are 
consistent with a flattening longwar d of 1700 A and a UV upturn 
shortward of 1700 A, as discussed bv lMaioIino et al.l II2004I) . There 
is no conspicuous kink in cither curves but they are clearly in- 
consistent with an SMC-like relative extinction curve. Note that 
the relative GIO or lMaiolino et al.l 1(2003) curves cannot be plotted 
because they are not defined above 3200 A. 

strumental problems may play a role for the NICS instru- 
ment which is a very low resolution instrument, making 
sky-line subtraetion complicated in the NIR. Conversely, 
NICS is a prism instrument and so is not susceptible to 
possible systematic slitloss correction errors in the way 
that the ESI and OHS/CISCO slit spectra may be. 

We have to accept that the available data on SDSS 
J1048-I-4637 can only lead us to inconclusive results 
about its extinction curve. This could be resolved with 
well-calibrated, simultaneous optical-IR spectrophotom- 
etry in the future. Unfortunately, the object is unobserv- 
able with X-shooter/VLT due to its declination. 

While we confirm that the flat extinction curve is in- 
consistent with a typical SMC curve, we cannot endorse 
the view that there is evidence for a kinked 'SN'-type 
extinction curve in this QSO, leaving no spectroscopic 
evidence for such a curve in any object to date. GRB 
071025 is the only currently uncontes ted case of such 
an extinction curve (jPerlev et al.ll2010l) . However, GRB 
071025 has only a photometric redshift at z ~ 5 (which 
impacts on the wavelength scale of the extinction curve) 
and the detection of a kink relies on _ff-band photometry 
being accurate to a few percent. 

6. DISCUSSION 

6.1. Inferring hot dust masses 

In the process of estimating extinction in high-redshift 
QSOs observed by IRAC / Spitzer, we have devised a care- 
ful method of estimating their intrinsic spectral slopes. 
The main systematic uncertainty in this method comes 
from estimating the parameter rj which controls the cor- 
rection for the emission lines entering the IRAC channels. 
Accurately estimating the intrinsic spectral slopes is also 



important for inferring hot-dust masses of these QSOs. 
Poor estimates of the spectral slope could lead to biases 
in the inferred hot dust component as measured essen- 
tially in the IRAC 8.0 ^m and 24 /im data. In Figure [7] 
we the refore compare our spectral slopes to those used by 
IJiang et al. (2006. 20101 in their estimates of the mass of 
hot dust in the targets. It is evident that their spectral 
slopes are generally lower than ours, which leads to an 
underestimate of the hot dust signal. In particular, the 
hot dust mass of SDSS J1306+0356 should probably be 
revised upwards because of its significantly high er spec- 
tral slope than that found by Jiang et al.l (| 20061) . Con- 
versely, the hot dust mass of J0840-I-5624 may have to 
be revised downwards, although we note that this object 
was found to have spurious negative extinction, due to 
either variability or problematic photometric data (see 
Section [ 



6.2. Correlations with other properties 

We compare our inferred spectral slopes and extinc- 
tions with other properties of the QSOs which might 
correlate with these values. In Figure [T3] we plot the 
black hole masses, cold dust masses and hot dust masses 
for the sample, to the extent that the data are available. 
There are no obvious correlations between the inferred 
spectral slope or extinction with either the mass of the 
black hole or the mass of hot or cold dust. However, fo- 
cusing in particular on the targets with sub-mm or mm 
detections or upper limits. Figure [T51 shows that the sys- 
tems with large amounts of cold dust (the filled symbols 
represent systems with sub-mm detections) also gener- 
ally have large black hole masses and hot dust masses. 
Moreover, there is an indication that systems with sig- 
nificant extinction observed with SCUBA have sub-mm 
detections. This is encouraging and suggests that a sig- 
nificant part of the signal we find is indeed due to extinc- 
tion and not an artifact of the data analysis or a thermal 
accretion disk. 

6.3. The case of SDSS J1044~0125 

Judging from Figure [H SDSS J1044-0125 appears 
to constitute a somewhat peculiar case because it is a 
BAL QSO with large amounts of inferred cold and hot 
dust and a large black hole mass. Yet, we only find 
evidence for very little extinction. We note, however, 
that if we assume a typical dust-to-gas mass ratio, and 
a spherical dust distribution in the host galaxy of the 
QSO, with angular diameter < 5" (Wang et al. 2010), 
we would_ anticipate Ay > 0.16 mag for the QSO. We 
observe ^1450 ^ 0.14 ± 0.24 mag which corresponds to 
Ay ~ 0.06 ± 0.16 for an SMC-like extinction curve (see 
Figure [T]) . These values are consistent within the un- 
certainties. Moreover, the dust distribution in SDSS 
J1044— 0125 is likely not spherically distributed. For ex- 
ample, the galaxy may be fiattened into a more disk-like 
shape. The relatively narrow wi dth of the CO (6-5 ) line 
relative to the [C II] linewidth (jWang et al.ll2013[ ) indi- 
cates complex gas dynamics. 

6.4. Caveats 

We seem to consistently find positive signals for the 
extinction and there is some correlation with dust signals 
found using other methods. Nevertheless, as we have 
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Figure 13. Black hole mass, cold dust mass and hot dust mass 
versus spectral slope and extinction (violet stars). Filled sym- 
bols represent QSOs with sub-mm detections, arrows in the mid- 
dle row represent QSOs with sub-mm upper limits. The rela- 
tive extinction values and spectral slopes for SDSS J1030-I-0524, 
SDSS J1044-0125, and SDSS J1048-I-4637 are based on the re- 
fined analysis of their spectra presented in thi s paper. Th e black 
hole masses and hot dust masses are from IJian g et al.| II2006I . 
12010 1) and M. Vestergaard (private communication). The cold dust 
masses are from SCUBA 850 iim observa tions tiPrLdde_v et al., 200^ : 
IRobson et"al|[2003: iPriddev et al.ll200a) . Note that the two sys- 
tems with Ai45o > 0.5 mag which have been observed with SCUBA 
have detections of large masses of cold dust. 

stressed throughout, our approach rehes on a number 
of assumptions, which we summarize below, along with 
additional possible systematic errors. 

• Likely, the most important caveat is the unknown 
shape of the intrinsic spectrum, as we have dis- 
cussed at some length. The assumption of a power 
law is a convenient one, for lack of a better credible 
and operational model for the intrinsic spectrum of 
aQSO. 

• We have estimated the effect of a thermal spectrum 
and shown that the effect is limited and perpendic- 
ular to the effect of extinction in the extinction- 
spectral slope plane (Figure [6]). Moreover, we have 
attempted to take out the effect of emission lines by 
dividing out a template excess continuum spectrum 
based on SDSS QSOs (Figure E|). On the other 
hand, the fake accretion disk extinction curves are 
quite reminiscent of the inferred extinction curves 
and we cannot rule out a small effect of a thermal 
spectrum in some of the systems. 

• We have determined the intrinsic power law slope 
using the IRAC /Spitzer detections in the 3.6 /Ltm 



and 5.8 ^m bands, avoiding the 8.0 ^m band which 
might be contaminated by hot dust emission and 
the 4.5 fim band because of Ha contamination. We 
note however, that the 5.8 /im band could in prin- 
ciple also be affected by hot dust emission in large 
quantities, although our analysis indicates that the 
effect will be small and only lead to an underesti- 
mate of the inferred extinction. A dominant contri- 
bution to the uncertainty in the extinction curves 
comes from the photometric uncertainty in the 5.8 
/im flux. 

• Moreover, the correction for emission line con- 
tamination of the 3.6 fim band is not trivial and 
we have adopted an estimate of the rj ratio from 
intermediate-redshift QSOs. A higher redshift tem- 
plate for more luminous QSOs would be highly de- 
sirable to constrain this ratio and its dispersion 
better. Fortunately, we have shown that the de- 
pendence of the inferred extinction on this r/ ratio 
is weak. 

• The IRAC/Spitzer data used to fix the intrinsic 
power law were not obtained simultaneously with 
other broad-band data or the spectroscopic data. 
Therefore, variability is a concern. Fortunately, it 
appears that variability is not a major issue for 
these very luminous QSOs (because of their very 
large black hole masses). Moreover, any variability 
is stretched in time by a factor '^ 7 due to the high 
redshift of the targets. However, we cannot rule 
out variability in individual cases. 

• Simple observational aspects such as thin cirrus 
during the spectroscopic observations, flux calibra- 
tion difficulties, or slit loss corrections might also 
lead to overall systematic errors of order 10-30%. 
Such errors could lead to very significant effects in 
the inferred extinction signal. 

• Foreground extinction in the Galaxy is fortunately 
small and has been corrected for. However, it is 
always possible that there may be extinction along 
the line of sight to an object, quite unrelated to 
the properties of the QSO. Foreground extinction 
along random line s of sight will be ab out an ob- 
served Ay ~ 0.1 (jMenardet al.ll2010f) . In addi- 
tion, there may be strong Mg II absorbers along 
the lines of sight to individual QSOs, such as that 
detected towards S PSS J1030-I-0524 at z = 2.78 
(jSimcoe et al.|[20Tl[) . In fact, we might expect an 
overdensity of Mg II absorbers towards very lumi- 
nous QSOs because of magnification due to gravita- 
tional Icnsing (as suggested for SDSS J1044-0125, 
see Shiova et al. 2002) . Foreground extinction may 
lead to extinction interpreted as host galaxy ex- 
tinction, and unusual extinction curves may arise 
from normal extinction curves contributing at the 
different redshifts, especially if the QSO extinction 
signal itself is weak. 

6.5. Theoretical SN extinction curves and the lack of 
evidence for a kink 

Despite the caveats discussed above, it does appear 
that there are significant extinction signals and that 
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these are related to the existence of dust in the QSO 
host galaxies. 

The overall shape of the obtained extinction curves pri- 
marily depend on the grain-size distribution, i.e., curves 
are typically flattened for grain-size distributions bi- 
ased towards large grains while specific features, such as 
for example the 1700-3000 A Todini-Ferrara-Maiolino 
bump, depends on the prevalence of specific grain species. 

The 'SN extinction curve' discussed in previous sec- 
tions is based on the dust composition and grain size 
distribution predicted by theoretical model s of du st for- 
mation in Type II SNe bv lTodini fc Ferraral |200l!) based 
on classical nucleation theorv (JFeder et all Il966l ') . The 
simplifications made include the assumption of uniform 
density and temperature as well as a full mixing of all el- 
ements in the SN ejecta. The predomina nt factors shap- 
ing t he Todini-Ferrara-Maiolino curve (jMaiolino et al.l 
120041 ) are small (« 10 A) silicate grains (primarily 
forsterite Mg2Si04), which cause the rise in the UV, and 
larger (si 300 A) amorphous carbon grains. The charac- 
teristic plateau at 1700-3000 A is ascribed to a minimum 
between two broad absorption features due to amor- 
phous carbon and a contribution from magnetite (F e3 04) 
(iMaiolino et al.l[2004 ). Bianchi fc Schneided (|2007f ) mod- 
ified the iTodini fc Ferrara (,20011) dust formation model 
and included a simple analysis of dust reprocessing due to 
a reverse shock penetrating the SN ejecta. This led to a 
shift of the grain size distribution to smaller grains (dom- 
inated by magnetite), however, maintaining the charac- 
teristic kink mentioned above, as well as the UV upturn 
below 1700 A. 

In the SN models for mixed and unmixed ejecta of 
iNozawa et al.l (|2003f ). who treat several physical process 
in more detail, non-carbon bearing dust is preferentially 
produced, with all components having grain-size distri- 
butions extendin g to l arger grains tha n predicted by 
iTodiiii fc Ferraral ()2001l ) . Consequently, iHirashita et al.l 
( 20051) predict rather featureless extinction curves for SN 
dust based on these models. Extinction curves from the 
mixed model show a steep rise in the UV whereas ex- 
tinction curves from the unmixed model are shallower. 
Using the same dust formation model but considering 
the change o f the grain-size dist ribution due to the re- 
verse shock (INozawa et al.ll2007[ ) flattens the UV slope 
([Hirashita et al.ll2008() as a consequence of efhcient de- 
struction of small grains. Such featureless, more shallow 
extinction curves w ere fa voured by GIQ, as w ell as by 
iKawara et aH (120111 ) and iShimizu et all (|201lD in z - 1 
ultraluminous infrared galaxies. 

A common drawback of predicted SN dust extinction 
curves is that they are precariously dependent on the SN 
dust formation models used (for a detailed review of SN 
dust models, see lGall et al.ll2011cl ). The models comprise 
different assumptions and simplifications about complex 
physical processes of the SN evolution and dust forma- 
tion leading to different results of dust grain properties. 
The size distribution is also very sensitive to the number 
of nucleation sites; more dust seeds will lead to smaller 
grains. Finally, the exti nctio n prop erties are usually de- 
rived by using standard iMid ()1908l ) theory for spherical 
grains. 

For exainple. iBianchi fc Schneideii (|2007D and 
iFallest et al.l (|2011[ ) have shown that the grain-size 
distribution in the dust-formation models sensitively 



depends on the assumed sticking probability, which is 
uncertain and commonly adopted to be unity. The 
sticking probability reflects the fraction of the col- 
liding condensible material which will stick together, 
subsequently building up macroscopic dust particles. 
Therefore, a larger sti cking probab ili ty res ults in larger 
grains. Moreover, IFallest et al.l (|201lD show that 
aspherical carbon grains tend to be larger. 

Finally, we note that none of the above predic- 
tions include the effects of the likely subsequent 
rapid reprocessing of dust in the inters t ellar medium 
of SN-rich galaxies (e.g. iJones fc NuthI f20Tlt iDraind 



200l iMichalowski et alj|2010al: IHirashita fc Kuol 120111: 
Hirashitall2012t lAsano et al.ll2013D ~ 

In this paper we have obtained the first non-parametric 
extinction curves for dust in high-redshift QSOs. No 
strong evidence for unusual extinctio n curves, such as 
that proposed bv IMaiolino et al.l (|2004( ). was found. This 
suggests that the overall dust properties at low and high 
redshift may be quite similar. We stress however that 
lack of a kink or other tell-tale extinction features do not 
preclude SNe from being dominant sources of dust in 
these systems. Given the strong sensitivity to detailed 
SN physics and unknown dust microphysics we should 
continue to seek model independent observational evi- 
dence for SN extinction. 
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